Chirped Biphotons and Their Compression in Optical Fibres 
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We show that broadband biphoton wavepackets produced via Spontaneous Parametric Down- 
Conversion (SPDC) in crystals with hnearly aperiodic pohng can be easily compressed in time using 
the effect of group- velocity dispersion in optical fibres. This result could foster important developments 
in quantum metrology and lithography. 
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One of the central problems in quantum optics is gen- 
eration of nonclassical light with given spectral and spa- 
tiotemporal properties. In particular, for the needs of 
quantum metrology and quantum lithography it is impor- 
tant to obtain two-photon wavepackets with small corre- 
lation times. Such wavepackets should naturally manifest 
a broad frequency spectrum. Several ideas have been put 
forward in this direction, all based on two-photon states 
produced via Spontaneous Parametric Down- Conversion 
(SPDC). Among them, one can mention prisms or diffrac- 
tion gratings introducing a frequency chirp [ij, SPDC 
in aperiodically poled crystals H, 0, SK and SPDC in 
crystals with temperature gradients pi]. However, a 
broad spectrum of two-photon light does not necessar- 
ily iniply small correlation times, although the inverse is 
true [y, 05 S|- This is similar to the fact that a broad- 
band pulse does not have to be short in time, although 
the spectrum of a short pulse is always broad. The spec- 
trum broadening introduced in Refs [H, 0, [H is in fact 
inhomogeneous; as a result, the two-photon spectral am- 
plitude in all these cases has a phase depending nonlin- 
early on the frequency. This phase (a frequency chirp [£]) 
makes two-photon wavepackets not Fourier transform- 
limited. Therefore, they are not short in time despite 
their broad frequency spectrum. As it was mentioned in 
Ref. , time compression of such two-photon wavepack- 
ets requires compensation for their frequency chirp. At 
the same time, the way to eliminate the chirp was not 
specified. 

In this paper we show that, under certain condi- 
tions, biphoton wavepackets can be made nearly Fourier 
transform-limited and hence compressed by injecting one 
of the photons of a pair in a standard optical fibre and 
exploiting the effect of group- velocity dispersion (GVD). 
No specially engineered fibres (for instance, with neg- 
ative GVD) are necessary. This suggests an easy way 
of achieving extremely short correlation times for two- 
photon light. 

Consider generation of two-photon light via sponta- 
neous parametric down-conversion (SPDC) from a cw 
pump in an aperiodically poled crystal. From the view- 
point of applications and for simplifying the calculation, 
it is convenient to assume that signal and idler photons 



are distinguishable, due to either frequency nondegener- 
ate or type-II phase matching. Below, we consider phase 
matching to be type-II, collinear, and frequency degener- 
ate, with idler (extraordinary) and signal (ordinary) ra- 
diations centered at frequency uoq. The two-photon state 
can be written as 



(1) 



where |c^)i(s) denotes the idler (signal) photon state 
with frequency uo. The two-photon spectral amplitude 
(TPSA) F(yt) determines all spectral and temporal prop- 
erties of two-photon light. In particular, its squared mod- 
ule gives the frequency spectra of signal and idler radi- 
ation, /s,i(c<;) oc \F{uj — a;o)P- Its Fourier transform can 
be called time two-photon amplitude (TTPA) [IQJ . 
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its squared modulegiving the second-order Glauber's cor- 
relation function [aQ: G^^\t) = |F(r)p. The TPSA is 
determined by the distribution of the quadratic nonlin- 
earity x(^) along the crystal [sl. fiol. [lH : 
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Here, L is the crystal length and ki^ks^kp are wavevec- 
tors of the idler, signal, and pump waves, respectively. 
Let the spatial dependence of the quadratic nonlinear- 
ity be x(^) = xoe*^^^^^^^^/^^ where the inverse grat- 
ing vector K has a linear dependence on the coordinate, 
K{z) = Ko - a{z + L/2) [3, 4, 9], and Kq provides 
quasiphasematching: ki{uo) + ks{uJo) — kp -\- Kq = at 
the center of the crystal. It is convenient to expand the 
wavevectors around the exact quasiphasematching fre- 
quency: 

1. 



ki = ki{uoo)-Kn^^k'/n^, 
ks = ks{coo)^k',n^^k'^n^. 



(4) 



Here, k^ ^ and k'/^ are the first and second derivatives of 
the dispersion law evaluated at cjq, related to the group 
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velocity and group velocity dispersion (GVD), respec- 
tively. 

Denoting D = k'^ — k[ and n = ^{k^/ + k^)^ we obtain 
the TPS A in the form 



iDn4 



, (5) 



where = z -\- L/2. 

Suppose that the spectrum is not too broad compared 
to the difference of group velocities of the signal and idler 
radiation, so that the condition 



D 



< 1 



(6) 



holds true. Then the wave vector mismatch can be writ- 
ten up to linear terms in frequency detuning and the 
TPSA becomes 

F{n) oc e-^^^* / ^ d^ xoe'''''^"'"^', (7) 

which yields, similarly to Ref. 
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In the case of large aperiodicity a, the spectral ampli- 
tude has a rectangular shape. Indeed, let us introduce 
the 'rectangle function' n(x, a, 6) = 1 for a < x < 6 and 
n(x, a, 6) =0 otherwise. Then, rewriting the integral in 
([7|) in terms of the rectangle function and applying the 
convolution theorem, we get 
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/ ax smciVx— \e 
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where sinc(x) = sm{x)/x. The first exponential term 
in the integral, exp{i ^^^ }, can be omitted if the typical 
scale of its variation is much larger than the sinc-function 
width, tt/DL. This is the case if the aperiodicity is large 
enough. 



a > 



(10) 



Note that condition (p!Q|) is well satisfied in Refs. [s*] and 
[1]. Then, the integral in Eq. ([9]) becomes 
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We see that the spectrum of SPDC in a crystal with 
linear K{z) dependence and large a is a rectangular func- 
tion of width Al] = The condition ([10|) means 



physically that the aperiodicity should induce a substan- 
tial spectrum broadening. Now we can explicitly write 
the condition ([6]) for the GVD of the nonlinear crystal to 
be negligible: 



nLa 



< 1. 



(12) 



For given a and L, this condition is realized if the sepa- 
ration of signal (ordinary) and idler extraordinary group 
velocities is large enough. 

Increasing the aperiodicity a, one can make the spec- 
trum as broad as desired. At the same time, this does 
not make the TTPA ([2]) short in time because, due to the 
nonlinear frequency-dependent phase factor in ([8]), the 
TPSA is not Fourier transform-limited. Here we would 
like to stress that because the squared module of TTPA is 
the second-order Glauber's correlation function, its width 
gives the correlation time of the biphoton, i.e. the bipho- 
ton entanglement time 0, 0] • It is this time that is impor- 
tant for two-photon effects such as two-photon absorp- 
tion, two-photon ionization, or up-conversion, and which 
can be measured for two-photon light using these tech- 
niques [l2[ . At the same time, coherence time of biphoton 
light is defined as the width of the first-order Glauber's 
correlation function, which is the Fourier transform of 
the spectrum \F(yt)\^. This is why coherence time does 
not depend on the phase factor in (|8|) and it is given by 
the inverse width of the spectrum |4| . 

The Fourier transform of TPSA ([7j) is easily obtained 
by introducing the rectangle function under the integral, 
extending the integration to infinite limits and using the 
convolution theorem. As a result, we get 
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with the amplitude being the same as in the case of a bulk 
or periodically poled crystal of length L. This means that 
the TTPA of an aperiodiocally poled crystal is as broad 
as in the absence of the aperiodicity a. 

Consider now propagation of the extraordinary pho- 
ton of the biphoton field through an optical fibre of 
length / with the inverse group velocity given by k'^ = 

£U=a;o and the GVD given by f^f = ^0U=a;o- Prop- 
agation through such a fibre leads to a phase factor 
eyip{i(k'j^Q + hifQ'^)l} in the two-photon spectral ampli- 
tude [6, 7]. The first term in the phase is linear in fre- 
quency and hence only shifts the two-photon wavepacket 
in time. The second term, being quadratic in frequency, 
can compensate for the TPSA chirp. This will happen 
under the condition 
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For a fibre with positive GVD, this condition can be sat- 
isfied for negative a, i.e., for the case where the poling 
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FIG. 1: (color online) Calculated spectrum of the signal/idler 
radiation (solid line). Dashed line: the corresponding spec- 
trum for a crystal with K = Kq, reduced in the vertical scale 
by 2 10^ 



FIG. 2: (color online) Calculated second-order correlation 
function of the biphoton in the case of aperiodic poling (black 
solid line) and periodic poling (green dashed line). 



period reduces along the pump propagation through the 
crystal. 

In the case of large aperiodicity a satisfying (p!Q|) . the 
resulting TTPA can be calculated analytically as the 
Fourier-transform of expression (pTj) with the quadratic 
phase term removed. Clearly, it has the form of a sinc- 
function with the width being almost aL^ times narrower 
than a periodically poled sample of the same length. 

Figure 1 shows the spectrum of SPDC radiation cal- 
culated for the case of aperiodically poled KTP crystal 
with L 0.8 cm, Kq = 2441.8 cm-^ and a = 1200 
cm~^, which corresponds (using the first-order quasi- 
phasematching) to the poling period varying from 18.47 
to 42.40 /im. The pump at 458 nm is y-polarized, as well 
as the idler radiation, and all three wavevectors, as well 
as the inverse grating vector, are directed along x. The 
dispersion dependencies are given by Sellmeier equations 
from Ref. [13], without any additional assumptions. One 
can see that the spectrum of the biphoton field is quite 
broad (from 800 to 1200 nm) and has nearly rectangu- 
lar shape. This is due to the fact that condition (p!Q|) is 
fulfilled very well. In fact, the aperiodicity leads to the 
spectral broadening of more than two orders of magni- 
tude. For comparison, the same figure shows the spectra 
of signal and idler radiation for a crystal with the same 
length but constant poling period K = Kq, correspond- 
ing to a = 0. Condition (p!2|) is reasonably satisfied, since 
1^1 « 0.16. 

Figure 2 shows the second-order correlation function 
calculated as the squared module of the Fourier trans- 
form of expression ([8]). For comparison, second-order cor- 
relation function of a crystal with constant poling period 
is plotted in the same graph. Clearly, both distributions 
have the same width, which means that the biphoton 
with the broadened spectrum has the same correlation 
time as the narrowband biphoton, i.e., it is not Fourier 
transform-limited. Evolution of the correlation time in 
an optical fibre is demonstrated in Fig. 3 for the cases of 
a > and a < 0, when only the idler photon is transmit- 
ted through the fibre. We see that at a > 0, propagation 
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FIG. 3: (color onhne) TTPA width at the output of the fibre 
versus the fibre length for a > (red Une) and a < (grey 
line). Circles denote the points for which the shapes of the 
second-order correlation function are shown in Fig. 4. 



through the fibre only broadens the biphoton wavepacket, 
while in the case a < 0, which is achieved by simply ex- 
changing the input and output faces of the crystal, the 
biphoton is compressed. For the calculation, we used 
the GVD of bulk fused silica jlj , because the waveguide 
contribution into GVD far from zero dispersion point is 
negligibly small [isl] . The value of tzf we used in the cal- 
culation was 1.359-10"^^^^. The largest compression of 
the biphoton wavepacket is achieved at the fibre length 
/ = 16.927 cm. The second-order correlation function in 
this case has a typical shape of squared sinc-function with 
the FWHM equal to 12 fs. At other lengths of the fibre, 
the two-photon wavepacket is broader. As the length of 
the fibre increases, the shape of G^^^(r) becomes similar 
to the shape of the spectrum, an effect that was studied 
in detail in Refs. [6l,[7[. Fig. 4 shows the shapes of G^'^\t) 
after the biphoton propagation through fibres of different 
length; the corresponding points are shown in Fig. 3. 

The model we have been using so far is based on the 
linear dependence of the wavevector mismatch on the fre- 
quency detuning from exact phasematching. This is valid 
for type-II or frequency non-degenerate type-I SPDC 
[lit , under the condition (p!2|) . Similarly, the dispersion 
law of the fibre in our consideration was described by a 
quadratic dependence, i.e., third-order GVD of the fibre 



4 



„4 

Z5 



-2-10123 
Delay x [ps] 




-3-2-1 1 2 3 
Delay x [ps] 




1.2| 

2 0.8: 



0.4 i 



-40 -20 20 40 
Delay x [fs] 


u 


-4-2 2 

Delay x [ps] 


4 




lilllttyiiilil 


_0.6 




f 






^ 0.4 




















0.2 










n 







-4 -2 2 4 
Delay x [ps] 



-4 -2 2 4 6 
Delay x [ps] 



FIG. 4: (color online) Second-order Glauber's correlation 
function of the biphoton with a < (a,c,e) and a > (b,d,f), 
after its propagation through a fibre of length 8 cm (a,b), 
16.927 cm (c,d), 50 cm (e,f) 
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Fibre length [cm] 

FIG. 5: (color onhne) TTPA width at the output of the fibre 
versus the fibre length for a > (red line) and a < (grey 
line) calculated without neglecting higher-order GVD terms. 



was neglected. In order to see the effect of higher-order 
GVD terms, we have performed exact numerical calcula- 
tion for the same case as considered above. The results 
(Fig. 5) show that the effect of compression is slightly re- 
duced but there still remains a significant narrowing of 
the TTPA, useful for applications. Even without any op- 
timization, the correlation time is reduced by more than 
an order of magnitude. An exhaustive study of this effect 
and a search for optimized parameters of both the crystal 
and the fibre will be presented in the nearest future [13] • 

An interesting feature of Fig. 5 is that at the output of 
the crystal, TTPA widths for the cases a > and a < 
are different. This can be explained as follows. Although 
at the center of the crystal signal and idler photons are 
generated with the same frequencies, biphotons gener- 
ated at the back face are non-degenerate. Due to GVD, 
there is a delay accumulated between the photons of a 
pair in the course of its propagation through the crystal. 



At a < 0, this delay is compensated by the one appearing 
due to birefringence and at a > 0, both delays add up. 

In conclusion, a biphoton whose spectrum is broadened 
due to a linear aperiodicity of the crystal poling can be 
compressed in time using normal GVD of an optical fi- 
bre. To describe the compression, it is sufficient to take 
into account first-order terms in the crystal dispersion 
dependence and second-order terms in the fibre disper- 
sion dependence. Exact calculation shows that higher- 
order terms reduce the compression but the effect is still 
present. It is worth mentioning an interesting result: the 
two-photon correlation time as well as its evolution due 
to the propagation of the biphoton through an optical 
fibre strongly depend on the position of the crystal. 

This work has been supported in part by the joint 
grant RFBR-Piedmont 07-02-91581-ASP, RFBR grant 
09-02-92003-NNS, MIUR (PRIN 2007FYETBY), Re- 
gione Piemonte (E14), "San Paolo foundation", NATO 
(CBP.NR.NRCL 983251), RFBR 08-02-00555a and the 
Russian Program for Scientific Schools Support, grant # 
NSh-796.2008.2. M. V. Ch. also acknowledges the sup- 
port of the Lagrange project of the CRT Foundation for 
her stay at INRIM. 



[1] A. Valencia et a/., Phys. Rev. Lett. 99, 243601 (2007). 

[2] S. Carrasco et al, Optics Letters 31, 253 (2004). 

[3] S. E. Harris, Phys. Rev. Lett. 98, 063602 (2007). 

[4] M. B. Nasr et a/., Phys. Rev. Lett. 100, 183601 (2008). 

[5] D. A. Kalashnikov, K. G. Katamadze, and S. P. Kulik, 

JETP, 89, 224 (2009). 
[6] M. V. Chekhova, JETP Lett., 75, 225-226 (2002). 
[7] A. Valencia et aL, Phys. Rev. Lett. 88, 183601 (2002); 

G. Brida et al Phys. Rev. Lett. 96, 143601 (2006) 
[8] D. Strekalov et a/., Phys. Rev. A 71, 041803(R) (2005). 
[9] Sometimes this dependence is called a linear chirp of the 

inverse grating vector; however, below we keep the term 

'chirp' for the frequency dependence of the phase. 
[10] A. V. Belinsky and D. N. Klyshko, Laser Physics 4, 663 

(1994). 

[11] G. Kh. Kitaeva, S. P. Kovalev, and K. A. Kuznetsov, 
International Journal of Quantum Information 7, 63 
(2009). 

[12] B. Dayan et al, Phys. Rev. Lett. 93, 023005 (2004); 

B. Dayan et al, Phys. Rev. Lett. 94, 043602 (2005). 
[13] V. Dmitriev et a/., Handbook of Nonlinear Optical 

Crystals (Springer Series in Optical Sciences, Vol 64) 

(Springer Verlag, Berlin, 1991). 
[14] H. Malitson, J. Opt. Soc. Am. 55, 1205-1208 (1965). 
[15] G. P. Agrawal, Nonlinear fiber optics, 3rd. ed. (Academic, 

San Diego, Calif., 2001). 
[16] Note that the same technique will be valid if the spectral 

broadening of two- photon light is caused by any physical 

mechanism leading to a linear variation of the wavevector 

mismatch along the crystal, in particular, by temperature 

gradients. 
[17] A. Meda et al.,m preparation. 



